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Abstract 

Taking the ground state rotational band in ^^Mg as an example, we investigate the impurity 
effect of A hyperon on collective excitations of atomic nuclei in the framework of non-relativistic 
energy density functional theory. To this end, we take into account correlations related to the 
restoration of broken symmetries and fluctuations of collective variables by solving the eigenvalue 
problem of a five-dimensional collective Hamiltonian for quadrupole vibrational and rotational 
degrees of freedom. The parameters of the collective Hamiltonian are determined with constrained 
mean-field calculations for triaxial shapes using the SGII Skyrme force. We compare the low-spin 
spectrum for ^'^Mg with the spectrum for the same nucleus inside ^^Mg. It is found that the A 
hyperon stretches the ground state band and reduces the B{E2 : 2^ — t- 0^) value by ~ 9%, mainly 
by softening the potential energy surface towards the spherical shape, even though the shrinkage 
eff'ect on the average proton radius is only ~ 0.5%. 

PACS numbers: 21.80.+a, 21.10.Rc, 21.60.Jz, 21.60.Ev, 23.20.Lv 
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I. INTRODUCTION 



Since the first discovery of A-liypernuclei by observing cosmic-rays in emulsion cfiam- 
bers [1], liypernuclei, wliicli are nuclei with one or more of the nucleons replaced with hyper- 
ons, have been used as a natural laboratory to study hyperon-nucleon and hyperon-hyperon 
interactions, properties of hadrons in nuclear environment, and in particular the impurity 
effect of hyperon in nuclear medium [2-4]. Due to the absence of Pauli's principle between 
the nucleon and the A particle, a A hyperon can probe deeply into the interior of nuclear 
medium and have important influences on its properties, including softening the equation 
of state [5], modifying the shape and size of finite nucleus [6], changing the nuclear binding 
and thus the driphnes of neutrons and protons [7] as well as the fission barrier heights in 
heavy nuclei [8]. 

In the past decade, many high-resolution 7-ray spectroscopy experiments using germa- 
nium detector arrays (Hyperball) have been carried out for A-hypernuclei [9] to understand 
the nature of A-nucleon interaction in nuclear medium and the impurity effect of a A on 
nuclear structure. In particular, the facilities built at J-PARC will provide an opportunity 
to perform hypernuclear 7-ray spectroscopy study with high precision by improving the 
quality of the secondary mesonic beam [10]. These facilities offer useful tools to study the 
low-lying states of hypernuclei, especially those of medium and heavy hypernuclei. To date, 
there are many experimental data not only on the single-A binding energy but also on the 
hypernuclear 7-ray spectroscopy that allow us to study the A-nucleon interaction, nuclear 
medium effects of baryons and impurity effects induced by a A hyperon in much greater 
detail [9]. 

The theoretical studies for the hypernuclear 7-ray spectroscopy are mainly performed 
with the cluster model [4, 11], few-body model [12, 13], and shell model [14]. The energy 
level scheme. Ml and E2 transition rates in low-lying states of light A-hypernuclei have been 
investigated with either a one-boson exchange potential or a parameterized spin- dependent 
A-nucleon interaction. Due to the numerical difficulty, the application of these models to 
medium and heavy hypernuclei is greatly limited. It is noted that, recently, the framework 
of few-body model has been extended to the case of five-body and used to study the energy 
levels of the double A-hypernucleus, ^a-^^ i^^]- 

The framework of nuclear energy- density functionals (EDF) is nowadays one of the most 
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important microscopic approaches for large-scale nuclear structure calculations in medium 
and heavy nuclei [16] and has already been extended to study hypernuclei [17-24]. Re- 
cently, both the non-relativistic Skyrme-Hartree-Fock (SHF) theory [25-27] and the rela- 
tivistic mean-field (RMF) theory [28] have been applied to study the impurity effect of A 
hyperon on the deformation of A-hypernuclei. The predicted energy surface is somewhat 
soft, in which case a large shape fluctuation effect of collective vibration might be expected. 
Furthermore, the static single-reference (SR) EDF is characterized by symmetry breaking 
(e.g., translational, rotational, particle number), and can provide only an approximate de- 
scription of bulk ground-state properties. Therefore, to calculate excitation spectra and 
electromagnetic transition rates in individual hypernuclei, it is necessary to extend the SR 
EDF framework to include collective correlations related to restoration of broken symmetries 
and to fluctuations of collective coordinates. 

In recent years several accurate and efficient methods and algorithms have been de- 
veloped that perform the restoration of rotational symmetries in 3D Euler space broken 
by the static nuclear mean field and take into account fiuctuations around the mean-field 
minimum [29-32]. The most effective approach to configuration mixing calculations is the 
generator coordinate method (GCM). Within these methods, the energy spectrum and elec- 
tromagnetic transition rates of low-lying excited states in both light and heavy nuclei have 
been successfully reproduced. However, these approaches are currently developed only for 
even-even nuclei, which cannot be extended straightforwardly to study the 7-ray spectra of 
single-A hypernuclei by simply adding hyperon degree of freedom. 

At present, the extension of 3D angular momentum projected GCM (3DAMP+GCM) 
method to single-A hypernuclei based on triaxial symmetry-breaking intrinsic states is still 
much complicated and its applications to medium-heavy and heavy nuclei would be computa- 
tionally demanding. As an alternative approach to the 5D quadrupole dynamics that restores 
rotational symmetry and allows for fiuctuations around the triaxial mean-field minima, a 
5D collective Bohr Hamiltonian (5DCH) has been formulated with deformation-dependent 
parameters determined by microscopic selfconsistent mean-field calculations [33-36] . In this 
work, we will construct a 5DCH with the parameters derived from the Skyrme-Hartree-Fock 
calculations for the nuclear core in a single A-hypernucleus and calculate the corresponding 
low-spin excitation spectra. The impurity effect of A hyperon on the collective motion of 
an atomic nucleus will be examined by studying the modifications of collective excitation 
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spectrum. In this way, we will in this paper concentrate on the modification of the core 
nucleus due to the addition of a A particle, leaving the evaluation of the spectrum of the 
whole hypernucleus ClS db future work. 

The paper is organized as follows. In Section II we present a brief outline of the 5DCH 
method and the Skyrme-Hartree-Fock approach for A-hypernucleus. The collective poten- 
tial energy surface, parameters in collective Hamiltonian as well as the resultant collective 
excitation spectra for ^^Mg and the same nucleus inside ^^Mg are given in Section III. A 
brief summary and an outlook for future studies are included in Section IV. 



II. THE METHOD 



A. Collective Hamiltonian in five dimension 



The collective Hamiltonian that describes the nuclear excitations of quadrupole vibra- 
tions, 3D rotations, and their couplings can be written in the form: 



H = T 



vib 



coll 



(1) 



where V^ou is the collective potential. The vibrational kinetic energy reads. 
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K = l 
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(3) 



with denoting the components of the angular momentum in the body-fixed frame of a 
nucleus. It is noted that the mass parameters -B/3/3, -B/37, -B77, as well as the moments of 
inertia X^, depend on the quadrupole deformation variables (3 and 7, 

= 45«/3^sin2(7-2K7r/3), k=1,2,3. (4) 

Two additional quantities that appear in the expression for the vibrational energy, that is, 
r = B1B2B3, and w = BppB^^ — B^^, determine the volume element in the collective space. 



The corresponding eigenvalue problem is solved by expansion of eigenf unctions in terms of 
a complete set of basis functions that depend on the deformation variables f3 and 7, and the 
Euler angles 0, 6 and ip [37]. 

The dynamics of the collective Hamiltonian is governed by seven collective quantities, 
that is, the collective potential Kou, three mass parameters -B/3/3, -B/37, and B^^, and three 
moments of inertia X^. These quantities are functions of the intrinsic deformations /3 and 7 
and will be determined by Skyrme-Hartree-Fock calculations with constraints on the mass 
quadrupole moments. 



B. Skyrme-Hartree-Fock approach for A-hypernucleus 

In Ref. [27], the computer code ev8 [38] of SHF+BCS approach has already been extended 
for the study of A hypernuclei. Therefore, in the following, we start from this approach to 
calculate the seven collective quantities in the 5DCH, as shown in Eq. (1). 

In the SHF+BCS approach for A hypernucleus, the total energy E can be written as the 
integration of three terms, 

E = J d^r[SN{r) + TA{r) + SNA{r)], (5) 

where £^Ar(r) is the standard nuclear part of energy functional, including both p/i-channel of 

the Skyrme force and pp-channel of the (5-force, as well as the kinetic energy density for the 

nucleons [38, 39]. Ta^t) = ta is the kinetic energy density of A hyperon. £^ArA(r) is the 

2mA 

interaction energy density between the A and nucleons given in terms of the A and nucleon 
densities [40], 

^NA = ^0(1 + ^^o)P^PN + ^{ti + t^){rApN + TnPa) 

+ lw^i'^PN ■ J A + VpA ■ J^)rArpA. (6) 

Here, pA, ta and Ja are respectively the particle density, the kinetic energy density, and the 
spin density of the A hyperon. These quantities are given in terms of the single-particle 
wave-function of A and occupation probabilities [39]. tQ,ti,t2,t^, and are the Skyrme 
parameters for the AN interaction. 
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The pairing correlation between the nucleons is taken into account in the BCS approxi- 
mation. The density-dependent 6-ioTce is adopted in the pp channel, 



I- po 

V{ri,T2) = -g — r— 



P(ri 



5(ri-r2), (7) 



Po 

where is the spin-exchange operator, and po = 0.16 fm~^. 

The HF equations for the nucleons and A are obtained by varying the HF energy (5) 
with respect to the corresponding single-particle wave functions and are solved by discretiz- 
ing individual single-particle wave functions on a three-dimensional Cartesian mesh. More 
details can be found in Ref. [38]. 

The method of quadratic constraints on the quadrupole moments of the nuclear density 
is used to find nuclear intrinsic wave functions (including the quasiparticle energies Ei, oc- 
cupation probabilities Vi, and single- nucleon wave functions ipi) corresponding to the desired 
quadrupole deformations [38, 41]. With these wave functions, one can calculate the moments 
of inertia in Eq. (4) using the Inglis-Belyaev formula [42, 43] 

^^ = T. ^''%~_,"e'^\ wj)\', (8) 

where k = 1, 2, 3 denotes the axis of rotation, and the summation of i,j runs over the proton 
and neutron quasiparticle states. 

The mass parameters Bfj^^{l3,'-f) are also calculated in the cranking approximation [44] 

with 

M(n),i.u{P. ^^ = 2^ [E-^EY ' '''''' ■ 

The mass parameters i?^,^ in Eq.(9) can be converted into the forms of B^^^B^^^B^^ with 
the following relations [45], 

B^^ = 5ooaoo cos^ 7 + 2502^02 cos 7 sin 7 -022022 sin^T, (Ha) 
B^^ = (^22022 - ^ooaoo) cos7sin7 5o2ao2(cos^7 - sin^7), (Hb) 
B^^ = i?22«22COS^7 - 2i?o2ao2COS7sin7 Soo«oosin^75 (He) 

where the coefficients 000,0025*222 are as follows, 

9^4^10/3 

O'OO = IQ^2 — ' '202 = 00o/v2, 022 = 000/2, (12) 
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FIG. 1: The collective potential VcoW, the rms radius of protons, the mass parameters -B/3/3, -B77, and 
the moment of inertia along the 1-axis Ii as functions of quadrupole deformation [3 for ■^^Mg and the 
nuclear core of ^^Mg from the Skyrme-Hartree-Fock+BCS calculations using the SGII force [46]. 
The AN interaction energy £na has been included (w) or excluded (w/o) in the collective potential 
FcoU for the nuclear core of ^^Mg. 

with ro = 1.2. 

The collective potential Kou in the collective Hamiltonian is obtained by subtracting the 
zero-point-energy (ZPE) from the total mean- field energy [35], 

Koii(/3,7) = i?tot(/3,7) - AKib(/3,7) - AKot(/3,7), (13) 
where Etot is the total energy for the nuclear core in A hypernucleus. We will investigate 
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two options, that is, those with (w) or without (w/o) the interaction part of energy Sna 
between the A and nucleons, 

i..o.4 (14) 

In the collective potential Vcou of Eq.(13), the vibrational ZPE, AKib is given by, 

AKib(/3,7) = ^Tr[>i(35A^(2)], (15) 

where A^(n),/i!/(/3, 7) is determined by Eq.(lO) with the mass quadrupole operators {fj,,^ = 
0,2) defined as , 

Q20 = - - and (^22 = - ■ (16) 
The rotational part of ZPE is a summation of three terms, 

AKot(A7)= ^^mm(/5'7), (17) 

with 

AV^u{f3,j) = -— — — . (18) 

where M.{n),fj.u{/3,l) is determined by Eq.(lO) with the intrinsic components of quadrupole 
operator defined as. 



Q 



2n 



III. RESULTS AND DISCUSSION 



—2iyz, /i = 1 

-2xz, /i = -1 (19) 
2ixy, /i = — 2 



Following Ref. [27], in the p/i-channel, we adopt the SGII parameterized Skyrme force [46] 
for the NN interaction, and the No.l set in Ref. [47] for the AN interaction. In the pp-channel 
for nucleons, we follow Ref. [48] to use g = 1000 MeV fm^ for both protons and neutrons. 
A smooth pairing energy cutoff of 5 MeV around the Fermi level is used. In the mean- 
field calculations, the mass quadrupole moments are constrained to the mesh-points in /3-7 
plane with /3 = 0, 0.05, 0.10, . . . , 1.20 and 7 = 0°, 6°, 12°, . . . , 60°. The A particle occupies 
the lowest single-particle state throughout the constraint calculations. In the following, 



8 



Mg A IVIg{core, w/o) ^ Mg(core, w) 




FIG. 2: The probability distribution pja{(3,^) in /3-7 plane for the 0^ state in ^^Mg (left panel), 
25 Mg (core, w/o) (middle panel) and ^^Mg (core, w) (right panel). 

we take ^^Mg as an example, and study the impurity effect of A hyperon by examining the 
changes of collective parameters and the resultant collective excitation spectrum and related 
obsevables. 

Figure 1 displays the collective potential V^ou, the rms radius of protons, the mass pa- 
rameters Bpi^, Bjj, and the moment of inertia along the 1-axis Ji as functions of quadrupole 
deformation /3 for ^^Mg and the nuclear core of ^^Mg. The AN interaction energy Sfqp^ in 
Eq.(5) has been included (w) or excluded {w/o) in the collective potential Vcow for the nu- 
clear core of ^^Mg [cf. Eq.(13)]. It is found that the A hyperon has negligible influences 
on the moments of inertia and mass parameters of the nuclear core. However, it can lower 
down the barrier in the neighborhood of spherical shape and make the energy curve stiffer 
at large deformed region. In other words, the A will reduce the collectivity of ^^Mg, where 
the AN interaction energy plays a major role. 

In Fig. 2, we plot the probability distribution p/^ in /3-7 plane for the 0^ state in ^^Mg 
and the nuclear core of ^^Mg, where the p/^ is defined as [49], 

P/. 7) = E I Kk^P^i) I I sin 87 1 , (20) 

K 

which follows the normalization condition, 

/ rf7P7„(/3,7) = l. (21) 

Jo Jo 



9 




FIG. 3: The difference in the probabihty distribution piaiPjj) of 0^ states between ^^Mg (core, 
w/o) and ^"^Mg (left panel) as well as between ^^Mg (core, w) and ^^Mg (right panel). 



E 3.10 



5 3.05 

o 

a 

0) 



01 

> 3.00 
n 



: o; 




5DCH(SGII) J 










(core, w/o) 


□ 

(core, w) 


Mg 


,Mg 


25.. 



- 101 :t 



0) 

□) 

100 iS 
c 

0) 

o 

0) 
Q. 



- 99 



FIG. 4: The rms proton radius for Mg and the nuclear core of ^ Mg. 

Here, \l/^j^(/3,7) is the collective wave function that corresponds to the solution of 5DCH 
in Eq.(l) and a = 1,2, ■ ■ ■ , labels collective eigenstates for a given angular momentum /. 
It is shown in Fig. 2 that the A shifts slightly the probability distribution of the 0^ state 
to the smaller deformation region. This effect can be seen more clearly from the changes in 
PiaiPil) for the 0^ state after the introducing of A hyperon, as shown in Fig. 3, where the 
differences in the probability distribution po.i(/3,7) for the nuclear core of ^^Mg and ^^Mg 
are plotted. Quantitatively, the average values of /3(7) are 0.54(20.0°) for ^^Mg and these 
values become 0.53(20.7°) for ^^Mg (core, w/o), and 0.52(20.8°) for ^^Mg (core, w). 

Moreover, it is also shown in Fig. 1 that the rms radius of protons is reduced by the A, 
in particular in the neighborhood of spherical shape. However, this shrinkage effect on the 
proton radius of ^^Mg is only ~ 0.5%, as illustrated in Fig. 4, where the rms proton radius 
in /3-7 plane for both the ^^Mg and the nuclear core of ^^Mg from the 5DCH calculations 
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FIG. 5: The low-spin spectra of the ground state band for the ^^Mg (b) and the nuclear core 
of ^^Mg (c, d) obtained by the five-dimensional collective Hamiltonian (5DCH) with the param- 
eters determined by the Skyrme-Hartree-Fock+BCS calculations using the SGII force [46]. The 
B{E2) values are in units of e^ fm^. The spectrum of ^^Mg is compared with the corresponding 
experimental data (a), taken from Ref. [50]. 

are plotted. 

Figure 5 displays the low-spin spectra of ground state band for the ^^Mg and the nuclear 
core of A^Mg. It is noted that the A stretches the spectra of ground state band. Comparing 
with columns (b) and (d), one finds that the A increases the excitation energy of 2^ state 
by ~ 7%. Moreover, it reduces the E2 transition strength B{E2 : 2^ Of) by ~ 9%, which 
is a little smaller than the values, 19(4)% or 16(6)% in ^Li [6]. 

IV. SUMMARY AND OUTLOOK 

The impurity effect of A hyperon in ^^Mg has been quantitatively studied in the frame- 
work of non-relativistic energy density functional theory that has been extended to include 
correlations related to the restoration of rotational symmetries and fluctuations of collec- 
tive variables by solving the eigenvalue problem of a 5DCH for quadrupole vibrational and 
rotational degrees of freedom, with parameters determined by constrained self-consistent 
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nonrelativistic mean-field calculations for triaxial shapes using the SGII Skyrme force. The 
low-spin spectra for ^^Mg in both free space and with the additional A have been calculated. 
It has been found that the A hyperon shifts the collective wave function of ground state 
to a smaller deformation region by softening the nuclear collective potential surface in the 
neighborhood of spherical shape. As the consequence of this effect, the spectra of ground 
state band becomes stretched and the excitation energy of 2f state is increased by ~ 7%. 
Moreover, the B{E2 : 2f — ?■ 0^) value is reduced by ~ 9%. However, the shrinkage effect 
on the average proton radius is found to be only ~ 0.5%. 

As pointed out in Refs. [26-28], the influence of the addition of A particle might be 
stronger in the relativistic mean-field approach. Therefore, it would be very interesting to 
extend this work to the relativistic case. In addition, to calculate directly the 7-spectra of 
single-A hypernucleus, one has to extend the current EDF based 3DAMP+GCM or 5DCH 
models for the odd-mass or odd-odd nucleus. Working along this direction is in progress. 
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Appendix A: Calculations of moments of inertia with the EV8 code 

In the ev8 code [38], the single-particle (s.p.) wave-function of /c-state $fc(r), discretized 
on a three-dimensional Cartesian mesh, is written in the 4-component form, 




(Al) 
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TABLE I: Parities of four components [a = 1,2,3,4) in single-particle wave function <&fc of 
/c-state with respect to the planes x = 0,y = 0, z = 0. The parity of /c-state is denoted as pk- 





X y z 




+ + Pk 

- - Pk 

- + -Pk 
+ - -Pk 



where {a = 1, 2, 3, 4) are real functions corresponding to the real and imaginary, spin-up 
and spin-down parts of The time- reversed state of $fc are determined by 



-v^« + ^^^^ 



(4) 



= r<i>A 



(A2) 



Therefore, the components in <l>^ are connected with the components in $fe by the following 
relations, 



(1) 



^(4) ^(3) 



(2) 



(A3) 



Together with the parities of the four components \E'[°'' in single-particle wave function 
$fc of A;-state with respect to the planes x = 0,y = 0, z = 0, as shown in Table I, the 
moments of inertia /i.2,3 in Eq.(8) can be simplified as 



i,j>0 



Ei + Ej 



i,j>0 



{uiVj - v.UjY r\j\-\'^ 



(A4a) 
(A4b) 



where the s.p. states i, j have the same parities {pi = pj) and the non-zero matrix elements 
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are determined by, 



(^I^iIj) 




+00 



dxdydz 



(2)^ 



(2)^i,T;(2) 



,{3)^ 




dxdydz 



(4)/ 




00 



dxdydz 



(3) 



(2)^ 



,1 

1 



(A5a) 



(A5b) 



(A5c) 



In the above equations, (k = x, y, 2;) denote the components of orbital angular momentum 
operator. 
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